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Abstract The homogeneous Boltzmann equation for 
inelastic Maxwell mixtures is considered to study the 
dynamics of tracer particles or impurities (solvent) im- 
mersed in a uniform granular gas (solute) . The analysis 
is based on exact results derived for a granular binary 
mixture in the homogeneous cooling state (HCS) that 
apply for arbitrary values of the parameters of the mix- 
ture (particle masses m,-, mole fractions c,, and coeffi- 
cients of restitution ay). In the tracer limit (ci — > 0), 
it is shown that the HCS supports two distinct phases 
that are evidenced by the corresponding value of E\ /E, 
the relative contribution of the tracer species to the to- 
tal energy. Defining the mass ratio fi = mi/ma, there 
indeed exist two critical values /Ujjfjg and ^^cg (which 
depend on the coefficients of restitution), such that 
Ei/E = for /LtgQg < fji < Mhcs (disordered or normal 
phase), while E\jE ^ for \x < fJ^^s ana Vor /Lt > /j^qq 
(ordered phase). 
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Granular assemblies depart from molecular systems 
not only from the difference of the length scales in- 
volved, but more importantly in that the interactions 
among constituents are dissipative PQ. In conjunction 
with the use of powerful experimental and numerical 
techniques, the application of non-equilibrium statisti- 
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cal mechanics to the field has yielded much progress in 
the last 20 years, whereas the questions were hitherto 
more centred on civil and mechanical engineering issues. 
Among the factors that explain this upsurge of interest 
for fundamental approaches, the pioneering work by I. 
Goldhirsch and G. Zanetti [2J, pertaining to the clus- 
tering instability in an unforced granular gas, should 
certainly be recognised. Isaac Goldhirsch subsequently 
became a leading figure in the granular matter commu- 
nity, much contributing to the improvement of meth- 
ods and understanding of models. As a tribute to his 
achievements and insights, we address here the partic- 
ularly simple situation of a mixture of grains, where 
spatial homogeneity is enforced, thereby discarding in- 
stabilities in the vein of the clustering phenomenon, but 
where non trivial out-of-equilibrium phase transitions 
take place. 

We consider a binary mixture of inelastic Maxwell 
gases at low density in the homogeneous cooling state 
(HCS). The corresponding set of coupled Boltzmann 
equations for the velocity distributions fi(y,t)(i = 1,2) 
then read 



d_ 

dt 



fi = J iA v \fi> fj]> 



(1) 



where the Boltzmann collision operator J%j[fi, fj] for 
dissipative Maxwell mixtures is 



Jn M/i,/j] = 

-Mvi)/i(v a )] 



da [a-V^vD/.K) 



(2) 



Here, ni is the number density of species i, fid is the 
total solid angle in d dimensions, and Qfjj < 1 denotes 
the (constant) coefficient of restitution for collisions be- 
tween particles of species i with j. Moreover, v' x = vi — 
fi 2 i (1 + a^) (CT-g)CT, V 2 = v 2 + Aii2 (1 + a^i) (S-g)?, 
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where g = vi — V2, <x is a unit vector directed along 
the centers of the two colliding spheres, and /iy = 
mi/(mi + rrij). 

The effective collision frequencies for collisions i- 
j are independent of the relative velocities of the collid- 
ing particles but can depend on space and time through 
its dependence on densities rii and granular tempera- 
ture T (see e.g. Ref. [3] for a discussion of this kinetic 
notion). They can be also seen as free parameters of 
the model. In previous works on multicomponent gran- 
ular systems [HEJEIj w rs was chosen to guarantee that 
the cooling rate for inelastic Maxwell models (IMM) be 
the same as that of inelastic hard spheres (IHS). With 
this choice, the collision rates u>ij are (intricate) func- 
tions of the temperature ratio T1/T2, which precludes 
analytical progress. Here, since our problem involves a 
delicate tracer limit, we aim at the simplest possible 
approach. Specifically, we assume that is indepen- 
dent of the partial temperatures Tj of each species but 
depend on the global temperature T = c\T\ + C2T2, 
Cj = rii/(ni + 712) being the mole fraction of species i. 
Thus, one considers the simple "plain vanilla" Maxwell 
model defined as w tJ = vcj, where v = AnVT is an 
effective collision frequency and the value of the con- 
stant A is irrelevant for our purposes. The form of u>ij 
is closer to the original model of Maxwell molecules for 
ordinary gas mixtures [7]. The plain vanilla Maxwell 
model has been previously considered by several au- 
thors [SlI MTUirn] in homogeneous problems pertaining 
to granular mixtures. 

In the absence of any external energy input, the 
granular temperature T{t) monotonically decays in time 
due to the inelastic nature of the collisions. We are 
here mainly interested in the time evolution of the par- 
tial temperatures Tj (or equivalently, the partial pres- 
sures pi — n,Tj). In the hydrodynamic regime (for times 
much longer than the effective mean free time v ), it 
is expected that all the time dependence of Pi is only 
through its dependence on the global temperature T(t) 
P] . The time evolution of T(t) is simply d t T = -(T 
where 



dnT 



m l I dvv 2 J ij {f l J,j 



(3) 



using v ~ vf and gets the usual Haff 's law [T3] 
T(0) 



T{t) 



[1 + K(0)*] 



(4) 



The partial pressures Pi can be determined by mul- 
tiplying both sides of Eq. (fTJ by rriiV 2 and integrating 
over velocity. Taking into account previous results [1] 
derived for dissipative Maxwell mixtures, one obtains 



0. 



where V is the column matrix 
V=( P \ 

\P*2 

and C is the square matrix 
'X + A u A 12 



£ = 



A 2 i A + A 22 



(5) 



(6) 



(7) 



Here, p* = Pi/nT and we have introduced the dimcn- 
sionlcss quantities 



2lv 



An = -^-(l-a? 1 )+^M2i(l+ai 2 ) 



1 - — (1 + CV12) 



4is 



w 12 Pi ,,2 f-, , „ \2 

a P2 



(8) 



(9) 



where w*j = Uij/v and pi = nivrii is the mass density 
of species i. The coefficients A22 and ^.21 can be easily 
obtained from Eqs. (J5J) and © by change of indices 1 <-> 
2. In addition, as pointed out earlier, the temperature 
T{t) behaves for long times as T(t) = T(0)e Ar where A 
is a nonlinear function of «y and the parameters of the 
mixture. 

After a certain kinetic regime lasting a few collision 
times, one expects that the reduced partial pressures 
p* and p* 2 = 1—Pi reach well-defined steady values p* s 
and P2 B , respectively. These steady values are obtained 
by solving the homogeneous equation CP = 0. This 
equation has a nontrivial solution if det C = 0. This is 
a second-degree polynomial equation whose largest root 
Amax = max(Ai,A2) = — (* governs the time evolution 
of the temperature in the long-time limit. Here, Ai and 
A2 are the solutions of the equation det C = 0, 



is the total cooling rate. In order to solve the temper- 
ature equation, it is convenient to change to a new 
time variable defined as r = J Q u(T(t'))dt' yielding 
T(t) = T(0)exp(-C*r) where C* = C/v. To find the 
relation between the "internal" clock (related to the 
average number of collisions suffered per particle) and 
the "external" time t, one integrates the relation for dr 



Ai 



-{An + A 22 ) T \/{An - A 2 2) 2 +4Ai 2 A 2 i 



2d 



The steady solution pi s is given by 
A12 



PlsW 



A12 - An - A' 



(10) 



(11) 
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Fig. 1 Plot of the temperature ratio T1/T2 versus the mass ratio 
m\/m2 for d = 3 in the case 022 = 012 = 0.8. The solid line is 
the result for IMM (Inelastic Maxwell Model) while the dashed 
line is for IHS (Inelastic Hard Spheres). 



They are given by 



Mhcs — 



a\2 



1+oS 



Mhos ~ 



"12 



l+c& 



(15) 

On the other hand, if the mass ratio /1 is smaller (resp. 
larger) than jUjjc S (resp. A'hcs)' tnen > ^2°' > - F° r 
elastic collisions (022 = ot\ 2 = 1), Mhcs = ^ an< ^ Mhcs — * 
00 and so A 2 is always larger than • 

The above results show clearly that there are two 
different regimes of behavior. When A^ ' > A^ , the 
tracer temperature T±(t) is enslaved to the granular 
temperature T 2 (t) and so, the temperature ratio asymp- 
totically reaches the steady state value p\ s (A 2 °' l )/ci. 
More explicitly, 



where A = A max . Consequently, for long times, the time 
dependence of the partial pressure p*(t) (or cquiva- 
lcntly, the energy ratio E\/E = £>*) can be written as 



Pi® 



A + Be-^-^) 1 - 



— (A2— Ai)vt 



where A and B arc constants depending on the initial 
conditions and the function p* s (A) is defined in Eq. 
pip . In conclusion, after a relaxation time of the order 
of KA2 — Ai)i^| _1 , the energy ratio p\ reaches the steady 
state value p\ s (A max ). As long as the mole fraction c\ ^ 
0, one has Ai 7^ A2 for any value of the mass ratio 
/i = mi/m 2 and the coefficients of restitution. It must 
be remarked that the results derived so far coincide with 
those previously obtained in the one-dimensional case 
®. 

Let us consider now the tracer limit (ci — > 0). In this 
limit, according to Eq. (|10[) . Ai and A2 become simply 



Ai -> A 



(0) 



-^ 21 (1 



«12) 



1 ^ 21 /1 1 1 

1 - — (1 +a 12 ) 



AO) 



a 



22 
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, (13) 



(14) 



The root gives the time behavior of the tempera- 
ture T\ of the tracer particles while \^ is associated 
with the time evolution of the granular temperature T 2 
of the excess component. For given values of the coeffi- 
cients of restitution 022 and a.\ 2 , it can be easily proved 
that A2 > ^1 if the mass ratio jj, lies in the range 
Mhos < M < Mhcs' wnere * ne critical mass ratios Mhcs 



and /^Hcg are obtained from the condition A^ 



A 



hm m 

t-»ooT a (t) 



2/Xl2A*2l(l + «12) 2 



4^2i(l + ai 2 ) [l 



"12)] - 1 



(16) 



On the other hand, if A^ ' > A^" 1 , then the combination 
(12) A12 — An — A j vanishes in the tracer limit so that, ac- 
cording to Eq. ((TT|) . the temperature ratio p* s (A^)/ci 
tends to infinity. This latter case corresponds to an ex- 
treme breakdown of the energy equipartition since the 
tracer particles are very energetic compared with the 
gas particles and the impurities essentially scatter off 
a static fluid background. The transition toward the 
heavy-impurity phase (i.e., when [i > Mhcs) was 
ready found by Ben-Nairn and Krapivksy [5] in their 
analysis on the velocity statistics of an impurity im- 
mersed in a uniform granular fluid. The light-impurity 
phase (i.e., when [i < Mhcs) ^ s n °l re P°rted in Ref. [8]; 



(0) 



it only appears when ai 2 > + a| 2 )/2. Thus, this 
new phase disappears (since Mhcs becomes negative) 
when a 12 = 0122 or when a% 2 < l/y/2. It must also be 
remarked that a similar non-equilibrium transition has 
been found for IHS [14] showing that, in the anomalous 
or ordered phase, the ratio of the mean square veloci- 
ties for the impurity and fluid particles Tim 2 /T 2 mi is 
finite (and so, the temperature ratio is infinite) even for 
extremely large mass ratios {m\Jm 2 — > 00). Although 
the transition to the heavy-impurity phase detected for 
IMM occurs in general for large mass ratios (for in- 
stance, Mhcs — 18.05 for a 22 = a\ 2 = 0.8), the transi- 
tion phenomenon found in Ref. |14j is less pronounced 
for hard spheres interaction since at a practical point 
one needs to consider much bigger values of the mass 
ratio for IHS to find the above transition. 

The expression (| 16(1 for the temperature ratio de- 
rived in the tracer limit (ci —> 0) agrees with the one 




Fig. 2 Plot of the order parameter E\/E versus the mass ratio 
mi/7712 for ct\2 = 0.9 and a\2 = 0.6. The hatched regions indicate 
the ordered phases while the arrows correspond to the critical 

9.83. 



mass ratios /^q S 
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^ 21 ft 1 ^2 
Li 2 = — ^-Mi2(l + ai 2 ) . 



(18) 



(19) 



Equation (fTTj) holds for Ai and A2 which are still func- 
tions of ci . To first order in c\ , they can be written as 



Ai(ci) w A^+A^ci, 



A a (ci)w4 0) +A^C!, (20) 



(i) 



where A^ and A^ are given by Eqs. (fT5|) and ([33]), 
respectively. The expressions of X^ and A^ can be 
obtained from their forms (jTOj) and (|11|) for arbitrary 
ci. After some algebra, one gets 



A« = - 



Mil^l2( 1 + a 12) 



2rf 



d 2 (A 



(0) 

2 



(21) 



obtained in Ref. [8] from the Boltzmann-Lorentz equa- 
tion. It appears that in general Ti/T 2 ^ 1, although Eq. 
p^|) shows that energy equipartition occurs when the 
mass ratio is given by fi cq = (l + a 2 2~2a 2 2 )/(l — a^), 
provided 022 7^ 1. Moreover, when the particles of the 
gas collide elastically (022 = 1), then Ti/T 2 = (1 + 
c*i2)/[2 + (1 — a\2) (1^21 / 1^12)]- This expression coincides 
with the one obtained jT5] for IHS. Beyond this case, 
the dependence of T\ /T2 on the parameters of the sys- 
tem in the normal phase is different from that of hard 
spheres [T^]. Figured] shows the dependence of the tem- 
perature ratio on the mass ratio for 022 = 0^12 = 0.8. 
In this case, /%cg — —0.055 and /l4ics — 18.05 and so, 
there is only heavy-impurity phase. The /x-dependence 
of the temperature ratio for IHS (with the same di- 
ameter for the tracer and gas particles) is also shown 
for comparison. We observe that IMM capture well the 
trends of IHS, except of course close to the mass criti- 
cal value where T1/T2 grows very fast with mi/7712 for 
IMM. This growing is less dramatic for IHS. Likewise, 
the plain vanilla approach exaggerates the features of 
the more refined Maxwell model alluded to above @J[5] 
Ej; the trends evidenced, though, appear to be robust. 

We now explore the physical consequences of the 
existence of critical mass ratios in the tracer limit. In 
order to analyze this point, we consider for instance the 
energy ratio E\/E = p\ s . For c\ 7^ 0, the expression 
of Ei/E is given by Eq. pip where the explicit depen- 
dence of An and A12 on ci, [i and ceij are given by Eqs. 
© and ©, respectively. If c\ — > 0, Eq. (fTTj) becomes 



Ei 
E 

where 

,(0) 



A 



ci 



(A$-Att)ci-A<$-X> 



(17) 



A^' = -A^ = =/i 21 (l + a 12 ) 



1 ^ 21 1 1 1 V 

1 - — (1 + a 12 ) 



X 2 1} = -^^12(1 + "12) 



1 



M12 



M21M12 (! + "12) 



d 2 (X { ° ] -A 2 u; ) 



(1 + CK12) 



It must be noted that if A = A^ in Eq. (JT7J) , then 
[according to Eqs. dHJ) and |[TSJ)] A$ + A 2 0) ^ and so 
the energy ratio E\/E = when c\ — > as expected. 



(22) 



However, if A = A^ 0) in Eq. (T7]), A^'+X^' 1 = and so, 
Ei/E ^ 0. More specifically, by taking the tracer limit 
in Eq. (fl7|) when A = Ai one gets 



(o) , ,(0) 



lim — =J- 

ci^o E 



A 



(i) 
12 



A 



.4 



(i) 



(23) 



= a 2 2 - 1 + 4/* 2 i(l + au) [1 - ^(1 + ai 2 )] 
a| 2 -l + 2/x 2 i(l-al 2 ) 

Note that although A\{ and A^ depend on the coeffi- 
cient of restitution an, the energy ratio is independent 
on collisions among tracer particles themselves. This 
means that one could neglect the Boltzmann collision 
operator Jn[/i,/i] in the kinetic equation of the one- 
particle velocity distribution function fi (Boltzmann- 
Lorentz description). This is quite a natural assumption 
when one analyzes the tracer problem. In addition, it 
is also usual to assume that the presence of tracer par- 
ticles does not affect the state of the solvent (excess 
component) and so, collisions of type 2-1 can be ne- 
glected (closed Boltzmann equation for the gas). On 
the other hand, in the heavy or light impurity phase, 
one needs to consider the contributions coming from the 
Boltzmann operator J 2 i[/ 2 ,/i] to get Eq. (|2"3"]l . This is 
clearly shown in the Appendix [X] where the expression 
(|23|) is derived from an alternative route. 

In conclusion, when A^cs < M < Mhcs' ^ nc tem- 
perature ratio T1/T2 is finite and so, the energy ra- 
tio Ei/E = 0. On the other hand, if \i < f%Q S or 
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Fig. 3 Phase diagram in the dissipation/mass ratio plane, for 
a fixed value of tracer-fluid inelasticity 012 = 0.9. The vertical 
dotted line indicates the cut corresponding to Fig. [2] 



/i > MhcS' ^h e temperature ratio diverges to infinity 
and the energy ratio becomes finite. This change of 
behavior is similar to an ordering process where the 
impurity is enslaved to the host fluid (E\jE = 0), or 
carries a finite fraction of the total kinetic energy of 
the system (Ei/E ^ 0). The latter situation can be re- 
ferred to as the "ordered" phase (extreme breakdown 
of the energy equipartition) while the first one can be 
coined "disordered" phase. Figures [5] and [3] illustrate 
the transition phenomenon found here. Specifically, Fig. 
[2] shows the energy ratio Ei/E as a function of the 
mass ratio m\jmi for 0:22 = 0.6 and ct\i = 0.9. In 
this case, according to Eq. (fT5|) . one gets /x[jcs — 0-041 
and /t%CS i= 9.83. It is apparent that, for asymptoti- 
cally small or large mass ratios, the contribution of the 
impurities to the total energy can be even larger than 
that of the host gas. Figure [3] shows a phase diagram 
in the {022, mi/m^j-plane at ct\i = 0.9. The light im- 
purity ordered phase appears for values of the mass 
ratio /i < 0.113 provided that 022 < 0.787 while the 
heavy impurity ordered phase is present for mass ratios 
/i > 5.487 in the complete range of values of the coeffi- 
cient of restitution 0^2 ■ Finally, it must be noted that a 
similar non-equilibrium phase transition has been found 
for a sheared granular mixture |16j . 



A Energy ratio in the ordered phase 

In this Appendix we will obtain the expression (23} of the energy 
ratio Ei/E in the ordered phase from the condition = £2 (this 
is the condition to determine the temperature ratio in the HCS 
1121 ) when the collisions among tracer particles themselves are 
neglected. Here, Q refers to the partial cooling rate associated 
to species i. In the case of IMM, the cooling rates Q have been 
exactly obtained in Ref. [4], In the tracer limit (ci — y 0), if one 
neglects the effect of collisions I-f , one has ~ £12 and £2 = 
C21 + C22 where 



C22 



C21 



l(0), 



M21M12 



1 _ a 22 
—V, 

2d 



(25) 



(26) 



where use has been made of the identity p% s = 1 — Pi B - Moreover, 
we are considering the ordered phase and so p\ is finite (and 
Ti/T is infinite). If one neglects the tracer collisions, the condition 
to get pj reduces to £12 = C22 + C21 Substitution of Eqs. (1211 ). 
1251 . and I I26H into the previous relation yields 



O ( A l - A 2 ) = ^21^12(1 + "12) -T" 

z z 1 Pl,s 

The solution to Eq. J27]l is 



Pi, 



a% 2 -1 + 4^21(1 + 012) [1-^(1+0:12)] 
a| 2 -l + 2M2i(l-af 2 ) 



(27) 



(28) 



Ci 



-A< C V 



^A*2i(l + 012) [l - ^-(i + "12)] 



(24) 



Upon writing this equation the explicit forms of A^ ' and 
have been considered. The expression II28H coincides with Eq. 
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